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g ■ Abstract 

fT) , The quark exchange model is a simple realization of an adiabatic approxima- 

tion to the strong-coupling limit of Quantum Chromodynamics (QCD): the 
quarks always coalesce into the lowest energy set of flux tubes. Nuclear mat- 
ter is thus modeled in terms of its quarks. We wish to study the correlations 
imposed by total wavefunction antisymmetry when color degrees of freedom 
are included. To begin with, we have considered one-dimensional matter with 
a SU (2) color internal degree of freedom only. We proceed by constructing 
a totally antisymmetric, color singlet Ansatz characterized by a variational 



$h ' parameter A (which describes the length scale over which two quarks in the 

system are clustered into hadrons) and by performing a variational Monte 
Carlo calculation of the energy to optimize A for a fixed density. We calculate 
the q — q correlation function as well, and discuss the qualitative differences 
between the system at low and high density. 
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I. INTRODUCTION 



With the advent of new, high-energy facilities for nuclear physics research (CEBAF, 
RHIC), the issue of identifying departures from a conventional, meson-baryon description 
of nuclear phenomena presses. One would naively expect that as the typical distance scales 
probed become short - relative to the radius of the proton - that the fundamental ingredients 
of the strong interaction, namely, the quarks and gluons, should become manifest in nuclear 
observables. This has not been the case. Indeed, unambiguous signatures of the formation 
of a "quark-gluon plasma" (QGP), hoped for at RHIC, have proven elusive. The fickleness 
of many newly-proclaimed signatures stem from the disparate pictures of observables in 
the quark/gluon or hadronic viewpoints. For example, J/ip suppression was heralded as a 
clear signature of QGP formation [Q], yet it seems that suppression can also be realized in a 
hadronic description of high density matter at high temperatures without plasma formation 
. The divorce of the two pictures is perhaps forced by the difficulty of treating quark 
confinement. We offer no new insight into this technical problem, but, rather, argue that it 
is useful to consider a simple model which interpolates between a hadron-based and quark- 
based description at low and high density. 

In this paper, we shall consider nuclear matter from the viewpoint of a constituent 
quark model. The particular model we use, the quark exchange model, has "natural" low 
and high-density limits: it behaves as a system of isolated hadrons at low density and 
as a non-relativistic Fermi gas of quarks at high density. Our purpose is to attempt to 
understand the structure of nuclear matter in terms of its quarks. Potentially, we would 
like to understand traditional issues, such as nuclear binding and saturation, from a quark 
model viewpoint. We anticipate that some "duality" may exist between the quark-gluon 
and hadronic descriptions. Indeed, the study of this quark-based description at nuclear 
matter density should help us to identify to what extent the descriptions are equivalent. 
However, this is not our sole purpose. In our, albeit simple, model, the structure of the 
system at fixed density is controlled entirely by the quark exchange dynamics, rather than 
by approximations to the hadronic and quark-based phases which are neither obviously 
consistent nor compatible. Consequently, we regard the study of the structure of the system 
with density as the major focus of our study. We hope to gain insight into the nature of a 
quark-hadron phase transition, as well as into the system's structure above and below the 
supposed transition point. In principle, our study also allows us to search for new collective 
excitations of the system, which could be novel in structure and unique to the quark-based 
description. However, we shall not explore such venues here. 

The model we use is the quark exchange model |4| . The quark exchange model is a non- 
relativistic, potential model, with many-body, confining forces. The many-body nature of 
the interactions allows the saturation of the hadron-hadron forces - there are no long-range 
hadron-hadron forces in nature - to be imposed by fiat. The original quark exchange model 
0j was constructed for the q 2 q 2 system in three spatial dimensions and in U(l) color. In 
this model, the U(l) limit is still confining, and the system studied shows rich behavior. For 
example, the symmetric spatial state possesses a weakly bound state. The binding energy's 
weakness, compared to the hadron confinement scale, is reminiscent of the deuteron. As a 
result, this model has been used in several previous quark model studies of nuclear matter 
and of finite many-quark systems [|5l-fr2]l. For the most part, the studies are confined to 
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one-spatial dimension and retain the restriction to U(l) color as well. These studies are 
of "spin zero fermions": spin is not explicitly treated, so that the spatial wavefunctions 
are required to be manifestly antisymmetric. This has the unfortunate effect of generating 
nuclear matter based on hadron clusters in parity-odd states. Horowitz and Piekarewicz 
have extended their work in one dimension || to three dimensions and have, in addition, 
considered an approximate extension to SU(3) color @. The studies, for the most part, 
rely heavily on variational Monte Carlo methods; however, the authors of Ref. M use a two- 
body reduction of the many-body Hamiltonian of Ref. |4| , so that they can use conventional 
many-body techniques. 

This work is concerned with the complete treatment of color degrees of freedom, and 
with the resulting structure of the system's ground state. In our model, the gluon degrees of 
freedom are passive - they act only to confine the quarks, so that they can not be excited. 
Only the quarks, then, carry an explicit color label; our treatment is "complete" in this 
context. As in previous papers, we shall continue to ignore the quarks' spin. As we do 
include color, though, nuclear matter is based on parity-even hadron clusters, as desired. A 
full treatment of color necessarily implies enumeration of all the exchange terms required 
for the construction of a fully antisymmetric state for a finite number of quarks. This 
necessity limits this initial study to SU(2) color and one spatial dimension, for tractability. 
We believe, however, that the techniques presented here can be extended to SU(3) color 
and three spatial dimensions. 

The simple model we use has many short- comings: it is not relativistic, even at high 
density, and there is no mechanism for the production and propagation of virtual qq pairs. 
Moreover, chiral symmetry has no meaning. Yet, we believe that the model study presented 
here gives useful insight into the nature of hadronic matter as a function of density. The 
model does include the effects of confinement and antisymmetry, and its virtue is that the 
structure of the system as a function of density is determined by this physics. To some 
extent, our results can be understood in terms of the interplay of these two effects. The 
qualitative insight we gain as a result should be relevant in understanding the predictions of 
richer models, as the effects we have included, confinement and antisymmetry, surely play 
no minor role. 

Section [H| summarizes the U(l) quark exchange model, and Sec. |l| describes its ex- 
tension to SU(2) color. We compare our model with other models of color in Sec. [TV[ , 
before proceeding, in Sec. |V|, to describe the Monte Carlo computations. Our results are 
summarized in Sec. [VI], and we conclude with an outlook towards future work in Sec. VT.I 



II. THE QUARK EXCHANGE MODEL 



Here we describe the quark exchange model. For two isolated quarks in one spatial 
dimension, we consider the Hamiltonian 

H = -\dl-\dl+ l -{x 1 -x 2 f . (2.1) 

The units have been chosen such that the quark mass m q and the oscillator frequency to 
are m q = h = uj = 1. We shall consider harmonic confinement exclusively, merely for 
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tractability. The potential term can be thought of as a flux tube which confines quarks 1 
and 2. For N quarks, we wish to generalize the potential term so that confining potentials 
"connect" the quarks to yield the lowest possible potential energy. We shall avoid any 
explicit discussion of color for the moment, and simply require that the quark "links" act in 
a pair- wise fashion. Following Horowitz et al. 0, the Hamiltonian for N quarks becomes 

1 N 

H=--Y,% + V{x 1 ...x N ), (2.2) 

i=i 

where the potential V{x\ . . . x^) is 

V(x\ . . . xn) = min Vpixi . . . xn) 
[p] 

such that 

N/2 

Vp(x) = v (x P 2i-i — x P 2i) (2.4) 
i=i 

with v(x P i — x P 2) = (x P i — x P 2) 2 /2. The P are elements of the permutation group Sjy, 
where P l is the 2 th member of the element P. In this way, all possible pair-wise confining 
potentials are sampled. 



(2.3) 



FIG. [I]. An illustration of quark exchange dynamics for 4 quarks. Two incoming hadrons 
with quark content (12)o and (34)o may exchange quarks - the strings connecting 1 to 2 and 
3 to 4 flip - to yield a (13) (24) outgoing state. 
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As the potential requires that the quarks be paired in a way that the global potential 
energy is minimized, the model is driven by quark exchange dynamics. That is, slight 
changes in the spatial configurations of the quarks may result in rearrangements of the links 
which confine them. The exchange of quarks between links is the only dynamical feature of 
the model. These dynamics are illustrated schematically in Fig. [1]. 

The advantage of the simple prescription of Eq. (|2.3|) is that cluster separability is guar- 
anteed. By "cluster separability," we mean that no long-range, residual forces exist between 
color-singlet clusters at large separations. Such color Van der Waals forces do not exist 
in nature [|Tj|, so that cluster separability is a desirable property. Yet in all models with 
pair-wise, confining forces, color Van der Waals forces (dipole-dipole interactions) do exist 



14j . The potential of Eq. (|2.3j ) succeeds in providing cluster separability because the forces 
are of an intrinsically many-body nature. 

In some sense, Eq. (|2.3| ) represents an adiabatic approximation to a strong coupling 
picture of quark-quark interactions. The gluon degrees of freedom, which are inert, act as 
links which connect the quarks to yield the lowest potential energy. We do not include the 
possibility of gluonic excitation in our model; the absence of low-lying qq states with exotic 
quantum numbers |15| supports this model choice. The interaction is adiabatic, as the links 
rearrange themselves instantaneously in adjustment to the configuration of the quarks. This 
picture relies on a separation of time-scale between the quark and gluon degrees of freedom 
in QCD. This is not unreasonable for heavy quarks; this is, fortunately, the same limit in 
which the non-relativistic quark model, which is also used here, has validity. For realistic, 
two-flavor quark matter, these approximations are not obviously robust; we adopt them out 
of necessity. 

The non-relativistic quark model has been justified a posteriori by its phenomenological 



success in reproducing a broad range of baryon properties |TB[. In addition, the U(l) q 2 q 
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model of Lenz et al. |4| evinces appealing qualitative behavior. As mentioned earlier, a weak 
bound state - compared to the confinement scale - exists in three dimensions. This scale is 
reminiscent of the deuteron. There also exist hidden channel resonances in the elastic phase 
shift, which occur at the hadronic inelastic channel thresholds. Yet, significant corrections 
to the adiabatic picture on which the Lenz et al. model is based likely exist. Existing studies 
focus on the N — N system, which suggest that the corrections are not small [|17],|]1|]. The 
simple Lenz et al. model does have interesting physics content, so that we shall proceed to 
explore it in an extension to quark matter in 577(2) color. 



III. MAGIC MODEL IN SU(2) COLOR 



Here we shall consider SU (2) color quark matter in one spatial dimension, although we 
shall continue to ignore the quarks' spins. We want to extend the model above in such a 
way that the ground state is guaranteed to be an overall color singlet and that the system 
reduces at low density to color-singlet hadron clusters. These constraints are not sufficient to 
specify the extension to 577(2) color uniquely. Starting with the £7(1) potential of Eq. ( p. 3D , 
one possible extension is 



V S U(2)(X 1 ...X N ) 



min Vp(x\ . . . Xn) 
[p] 



(a + 1 — aP[ 



(3.1) 
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where P[p min ] denotes a color singlet state projector for the minimum energy pairing and a 
is a large and positive constant. For the pairing 12 34 56 78 . . ., 

P [Pmm]o = |(12) (34) (56) (78)o---)((12) (34) (56)„ (78) • • ■ | • (3.2) 

The potential choice in Eq. ( |3.1[) yields the required low-density limit of color-singlet hadron 



clusters. The constant a must be large re 
of the potential, that is, large relative to 



ative to the maximum value of the spatial part 
min[p] Vp(xi . . .£jv) , in order for the prescrip- 



tion to be sensible. The simulations are performed for quarks in a box of length L under 
(anti)periodic boundary conditions, so that the spatial potential is bounded for fixed L. 
The state |(12) (34) (56)o (78)o • • •) is a manifest color singlet, though it is only a single 
component of the full antisymmetric wavefunction for the quarks. One must begin, then, 
with a Slater determinant for the initial wavefunction in order to guarantee antisymmetry 
under the exchange of any two quarks. This is no different from the U(l) calculation. The 



SU(2) prescription given by Eqs. (3.1) and (|3.2|) neglects any contribution to the potential 
from paired triplet states which are coupled to an overall singlet. The so-called "hidden 
color" forces are zero, then, in this prescription. 

Computing observables with the color projector of Eq. Q3.2p is cumbersome at best. 
In addition, the introduction of the constant a is unappealing. We shall choose, then, a 
different SU(2) extension for use in this paper. We shall retain the original U(l) form of 
the potential, Eq. ( |2.3| ), and impose the needed constraints on our variational Ansatz. The 
constraints we require are as follows. 

• The total wavefunction must be antisymmetric under quark exchange. 

• The total wavefunction must be a color singlet. 

• The total wavefunction must separate into color singlet hadron clusters as pairs of 
quarks are pulled apart. 

By imposing these constraints directly on the variational Ansatz, we have declared that the 
flux tube dynamics are such that these properties are simply guaranteed. For that reason, 
we call our extension to SU(2) C the "magic" model. This model choice implies that hidden 
color forces are now finite. The treatment of hidden color here is compared with previous 
treatments of color in the quark exchange model in section [TV]. 

We shall now proceed to construct the desired variational Ansatz, and shall illustrate the 
procedure to come by considering the following Ansatz for four quarks: 



(3.3) 



The notation |(12)o) denotes the pairing of quarks 1 and 2 into a color singlet. The functions 
/ FG are the spatial wavefunctions appropriate for a Fermi gas. Upon putting the quarks in 
a box of length L and choosing anti-periodic boundary conditions, we have 

A F 234 = cos T( X 1 + X2-X3- X A ) . (3.4) 



Atrial 0&1, X 2 ,X 3 ,X 4 ) = exp| 


—X(xx 


- x 2/ 


2 - A(x 3 


- *4) 2 ] 


fFG 
Jl234 


|(12)o(34) > 


— exp 


-\(X! 


-X3] 


2 - X(x 2 


- ^) 2 ] 


fFG 
il324 


|(13) Q (24)o> 


+ exp 


—X(xi 


- xa, 


2 - A(x 2 


- *3) 2 ] 


fFG 
J 1423 


|(14) Q (23)o> 
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In SU(2) C , two quarks can be put in a momentum state +k. The variational parameter A 
describes the clustering of the quarks into hadrons. Equation (|3.3|) has several important 
properties. 

• The wavefunction is explicitly antisymmetric under quark exchange. 

• The A = wavefunction is equivalent to a Fermi gas Slater determinant. 

• The wavefunction separates into color singlet hadron clusters in the limit of large 
hadron-hadron separations, that is, for example, in the \x\ — a^l, \%3 — — ► oo limit. 

The above wavefunction Ansatz can be readily generalized to iV quarks. One starts by 
generating all distinct sets of quark pairs - in the case of four quarks, there are three such 
sets (12), (34), (13), (24), and (14), (23). One then determines the "symmetry" of each 
pairing with respect to one given pairing. That is, one counts the number of permutations, 
m, required to bring the pairing considered to the original pairing and then (— l) m determines 
whether that pairing is even or odd. That sign is the over-all sign with which its contribution 
is added to the wavefunction. The contribution that a particular pairing (12) (34) (56) (78) ■ ■ • 
makes is exp[— X(xf 2 + ^34 + + x 78 )]f (si 2 , s 34 , s 56 , s 78 . . .), where = X{ — Xj and 
Sij = Xi + Xj. The function / FG is a symmetric function of its arguments, constructed from 
cosines of the N/4 momenta and all possible pairs of the s^'s. The allowed momenta in the 
box are odd multiples of ir/L, as we choose to enforce anti-periodic boundary conditions. 
Thus, the function / FG for the configuration 12345678 ... is 




Note that "per" denotes a permanent. A permanent is the symmetric version of a deter- 
minant, so that the various terms are just added together. In constructing the function 
/i2345678...! we sum over an possible sets of distinct pairs of the hadronic coordinates s^, 
without regard to their order. A set of hadronic coordinates is itself labelled u Ph" , and 
there are (A^/2)!/(2 iV / 4 (A^/4)!) different sets of distinct pairs for a given Ph. A member of 
Ph is denoted by P^. For a given P^, we construct the permanent of the allowed momenta 
and the s^- coordinates. In writing Eq. (|3.5|) , we have assumed that the number of quarks 
is a multiple of 4. It is only in this case that the plane wave states for each free quark 
can be collapsed to an explicitly real wavefunction. Including the Fermi gas wavefunction is 
cumbersome, but it is absolutely necessary. In its absence, the variational procedure is not 
well-posed, as the A = wavefunction is trivial. 

In this manner, the total wavefunction can be constructed, and it satisfies all the points 
highlighted above for four quarks. Specifically, for eight quarks, we have: 

tt£ai(zi, . . . , x 8 ) = exp[-X(x 2 l2 + x 2 u + x 2 56 + x 2 78 )} ff 2 % m78 |(12)o(34)o(56)o(78) > 

- exp[-A(o4 + x 2 24 + x 2 56 + x 2 8 )} f™ 45m |(13)o(24) (56)o(78)o) (3.6) 

+ exp[-A(x? 4 + x 2 23 + x 2 m + x 2 8 )} /f« 35678 |(14)o(23) (56)o(78)o) 
+ . . . (102 terms) , 



7 



where 

7T 37T 37T 7T 

A F 2345678 = COS ^ ( S 12 - s 34) COS -p(«56 ~ S 78 ) + COS -J~(Sl2 ~ S34) COS -(s 56 - S 78 ) 

7T 37T 37T 7T 

+ COS — (S12 - s 56 ) cos -p(«34 - S7s) + cos -y-(si2 - sse) cos — (s 34 - s 78 ) (3.7) 

7T 37T 3"7T 7T 

+ COS — (S12 - S 78 ) COS -^(^34 ~ S 56 ) + COS -^(«12 ~ S 78 ) COS -(s 34 - S 56 ) . 

The difficulty with this procedure is that the number of distinct pairs of quarks grows 
greatly as N increases. There are N\/(2 N ^ 2 (N/2)\) terms in the N quark wavefunction, so 
that the calculation quickly becomes unmanageable. There are 105 terms in the 8 quark 
wavefunction, but there are 10395 terms in the 12 quark wavefunction. 

Enumerating all the exchange terms by brute force is not practical, nor is it necessary. 
Our interest is not in the wavefunction itself, but, rather, in the observables associated with 
it. We are interested in the matrix elements (\l/|0|\l/), where O is a symmetric - and color- 
independent - operator. For example, O = T = — ; O can be the kinetic energy 
operator. It turns out that permutation symmetry greatly reduces the number of non-trivial 
terms. Consequently, for 8 quarks, say, the matrix element of the symmetric operator O can 
be written 

« 9 ial | 6 |< 9 ial > OC J{d Xl }f 12U5678 6 [(1)/ 123 45678 - (12)^/14325678 

+ (32)-/i4365278 + (12) -/l4325876 — (48) -/l4365872 , (3.8) 



where /12345678 = /12345678 ex P[~ -M 3 ^ + x 34 + x \§ + x 7s\- F° r the matrix elements of an 
operator O of the type above, one need only count occurrences of topologically distinct 
terms. These are terms which can not be brought into each other under rotation of the 
coordinate indices. For example, for eight quarks, starting with 12 34 56 78, there are 

"singles" «-> 14 32 56 78^ "[1]" 
"doubles" <-> 14 36 52 78^ "[2]" 
"triples" «-> 14 36 58 72^ "[3]" 
and 

"disconnected doubles" «-> 14 32 58 76^ "[1+1]" . 

The terms in brackets have been introduced as a convenient notation, so that we can extend 
this construction to iV quarks. We use "singles" or "[1]" to denote the single exchange of 
quarks between two clusters, whereas "doubles" or "[2]" denotes the exchange of quarks 
between three clusters such that one quark in the first cluster ends up in the third. This is 
in contrast to a "[1 + 1]" exchange, in which two single exchanges occur between two separate 
pairs of clusters. In Eq. (fTST), the numbers in parentheses are the number of terms of that 
type in the original wavefunction, so that their sum is 105. The fractions are the color matrix 
elements for terms of that type, note that ((12) (34) (56) • • • |(12) (34) (56) . . .) = 2 N ' 2 . 
Since we do compute the color matrix element, we need only consider permutations between 
one member of a particular pair of quarks; in the above, for example, we have kept the odd 
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quarks fixed and considered exchanges of the even ones. The sum of all the coefficients, 
then, including the color matrix element, is (AT/2)!, or 24 for the eight quark case. 

For iV quarks, one can write the general form of the coefficient, given its topological 
structure. Moreover, one can enumerate all the possible topological structures for a fixed 
number of quarks. In order to discuss this, we shall use the u [i + j + k + ■ - •]" notation 
defined above. 

The distinct topological terms for a fixed number of quarks are given by [i + j + k + ■ ■ ■], 
where the range of i,j,k,... satisfies the following bounds. That is, % ranges from 1 to 
the upper bound £ = N/2 — 1 and j, k, . . . range from to £, subject to the constraints 
that i + j + k + • • ■ (k terms) < £ — k + 1 and i < j ' < k < . . . when j,k, . . . ^ 0. 
For eight quarks, the terms are [1], [2], [3], [1 + 1]. For twelve quarks, the terms are 
[1], [2], [3], [4], [5], [1 + 1], [1+2], [1 + 3], [2 + 2], [1 + 1 + 1]. The formula for the overall 
coefficient of the [i + j + k + • • •] term can be expressed in closed form. That is, 



r - 1 



"I) 



i+j+k+- 



(i + l)(j+l)(k+l) 



(n H ) ■ ■ -{n H - i - j - k) , (3.9) 



where n# = N/2, and s is a symmetry factor, s is determined by the number of repeated 
exchanges. That is, if [i + j + k + • • •] has I terms of type i and J terms of type j, and all 
the other terms are distinct, then s = I\ J\. For example, with [1 + 1], s — 2. Thus, we have 
a strategy for the computation of the (^trSil^l^tSii) matrix element: 

(<L\6\<^ oc J{d Xl }f 12U ... N 0^ ck (3.10) 

where 

*trick = /l234-JV + C[i+3+k+~]f[i+j+k+-] ■ (3-H) 

i,j,k,... 

The notation J2' means that the sum is restricted. As above, we require that i < j < k . . . 
when j,k,...^0 and i + j + k + • • • (re terms) < N/2 — k. Note that i starts at one, whereas 
j, k, . . . begin at zero. In Sec. [V], we discuss the stochastic evaluation of Eq. ( |3.10| ) and 



the resulting computation of the total energy and the two-body correlation function with 
density. 

We can compute the (^trteil^l^triai) ma trix element in another manner. That is, 

<*fii|6|*fii> « J{dx t } «f al ) nn o (*£i) nu • (3.12) 

Here "|" and "J," denote the two color states of the model, and (^triaduU- denotes a 
projection of the full wavefunction on the (Tit! • • • I state. We have written the integrand in 



Eq. ( |3.12D as a manifestly positive definite quantity. This is highly convenient for a stochastic 



evaluation of the integral, though it is at the cost of a more complicated wavefunction. This 
wavefunction contains (N/2)\ terms, but it is simply related to \I/^ k . Starting with ^^k, 
we enumerate all the terms of a particular topological class. Thus, for 8 quarks, in place of 
the term with 
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[1] «-> 14 32 56 78 



(3.13) 



we explicitly include all 6 terms of that type. That is, 



14 32 56 78 
16 34 52 78 
18 34 56 72 
12 36 54 78 
12 38 56 74 
12 34 58 76 

In this manner, (^tr^TUi- can ^ e reacm y generated. 

We shall study the benefits and disadvantages offered by Eq. ( |3.1U| ) and Eq. ( 3.12| ) in 
Sec. [VJ 



IV. COMPARISON WITH OTHER MODELS 



Here we compare the model defined by Eq. ( |3.3|) with previous attempts to incorporate 
SU(N) color. To our knowledge, our model study is the first complete treatment of SU{2) C 
in quark matter. To make our treatment of hidden color forces clear, we first consider 
the original SU(N) C study of Lenz et al. [Q in the q 2 q 2 system. For a particular pairing 
(12) (34), for example, they postulate that confining forces of differing range may reside in 
the orthogonal sectors |(12) (34) ) and |(12) 1 (34) 1 ) . The sector in which the individual 
quarks are paired to triplets before being paired to an overall singlet is the "hidden color" 
sector. 

The color kets of Eq. ( |3.3j ) can be recast into the framework described above. That is, 
the kets |(13) (24) ) and |(14) (23) ) can be expressed as linear combinations of |(12) (34) ) 
and | (12)i(34)i)o- Following Lenz et al. [|J, we introduce the coordinates 



(4.1) 



These are useful reduced variables, as they describe the relative separation of hadrons with 
a particular quark content. With Eq. ( |4.1| ), we can rewrite Eq. ( |3.3| ) as 



X = 


1, 

2^1 + £3 


- x 2 


— X4 


y = 


I, 

2^1 + X 2 


- x 3 


— X4 


z = 


1 , 

-{Xi + x 4 


- x 3 


- x 2 



exp[-2A(x 2 + z 2 )] f? 2 % - -exp[-2A( 2/ 2 + z 2 )} 



lexp[-2A(x 2 + 2/ 2 )]A^3 



l(12)o(34) ) 



+ 



\ exp[-2A(y 2 + z 2 )] /J& + \ exp[-2A(a; 2 + y 2 )] f?g 3 



|(12)i(34)i) (4.2) 
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where |(12)i(34) 1 ) = | (12) 11)|(34) 1-1) + |(12) 1-1) | (34) 11) — | (12) 10) | (34) 10). Note that 
in |(12) IJ) / and J denote the z component of the color of quarks 1 and 2, respectively. We 
see manifestly, then, that hidden color states are equivalent to the coupling of the other pairs 
of quarks to singlets, and, thus, are trivial. Indeed, in the absence of active gluon degrees of 
freedom, this is generally so. In the channel with the pairing (12) (34), the relative separation 
of (12) and (34) is determined by y - so it is the "channel" variable. The exchange terms in 
the color singlet and hidden color sectors have the same behavior in y, so that the spatial 
range of the confining forces in the two sectors is the same. Lenz et al. [|j] considered the 
q 2 q 2 system as a function of a parameter specifying the relative spatial range of the color 
singlet (comprised of color singlets) and hidden color sectors. Our model of hidden color 
forces maps to one particular case in their study, although there is no especial physics reason 
for our choice. 



The full SU(2) treatment of Sec. |T| can also be compared with the approximate treat- 
ment of Horowitz and Piekarewicz 0. Their work was in SU(3) color and three spatial 
dimensions; we now consider this approach in SU(2) C for purposes of comparison. In their 
treatment of color, the quarks are given fixed colors. (Watson || has also used this starting 
point.) For this reason, in the discussions to come, we will refer to this approach as the 
"painted" model. In constrast to Eq. ( |2.3| ), the potential is no longer color-blind, but, rather, 
consists of the minimum energy pairing of quarks of differing color. The construction of a 
totally antisymmetric, color singlet variational Ansatz is straightforward. That is, 

*2i = exp[-AVk B ] v^ G < G , (4.3) 

where 

N/2 

[l,m\ i j=1 

such that v(x r Ll — x b M1 ) = (x r Ll — x b Ml ) 2 /2. The L and M are elements of the permutation 
group S(n/2), and L l is the i th member of the element L. We have presumed that there are 
equal numbers of red (R) and blue (B) quarks and that their coordinates are labelled by the 
appropriate superscripts. Note that \I/^ G is the Slater determinant for the red quarks. The 
following convenient form is equivalent to a Slater determinant: 



N 

7T 



*i G = n ^ 

n.m odd 



\X n X Ti 



(4.5) 



This method is simple to implement and is able to treat many quarks with ease. We 
will compare the full SU(2) C results to the simple approximation described here for two 
reasons. First, we wish to examine the extent to which the predictions of the painted model 
are robust. Its primary advantage is convenience; we can not expect the painted model 
to be exact as all color flipping transitions, i.e., RB <-> BR, are neglected. However, the 
differences might not be large. Second, we hope that any departures between the full SU (2) c 
and painted calculations will give us greater insight into the nature of the correlations in 
the full calculation. 
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We can begin to understand the difference between the painted (approximate SU(2) C ) 
and magic (full SU(2) C ) models through consideration of their variational Ansatze for four 
quarks. That is, for the magic model, 

(TITI |*m) = exp [-A[(ar 1 - x 2 ) 2 + (x 3 - x 4 ) 2 ]] cos f (xi + x 2 - x 3 - x 4 ) 

- exp [— X[(xi - x 4 ) 2 + (x 2 - x 3 ) 2 ]] cos j(xi + x 4 - x 2 - x 3 ) (4.6) 

whereas, in the painted model, 

(TITI I^p) = exp -Amin[V(i2),(34), V(i4),(23)] sin —(&! - x 3 ) sin ^(x 2 - x 4 ) (4.7) 

where quarks 1,3 are red ("up") and quarks 2,4 are blue ("down"). The Hamiltonians of 
the two models differ, but it is still useful to compare the structure of the Ansatze as two 
quarks of the same color approach each other. In this way, we can understand more about 
the dynamical content of the two models. Letting x\ — x 3 = e be small, we have 

(TITI l^tl) = exp [-A[(x 2 - x 3 f + (x 3 - x 4 ) 2 )] {-2ef sin f {x 2 - x 4 ) 

+ 2X(x 2 - x 4 )e cos f (x 2 - x A ) + 0{e 2 )} (4.8) 

for the full SU(2) C calculation and 



7T . 7T 



(TITI I^p ) = exp [-X[(x 2 - x 3 ) 2 + (x 3 - x 4 ) 2 ]\ {e- sin -(x 2 - x 4 ) + 0(e 2 )} (4.9) 

for the painted model. Equation ( |4.(j| ) contains an additional term in leading order in e. 
This term does not have the same manifest L dependence as the term which both Eqs. ( |4.8|) 
and (|4.9| ) share. Indeed, this additional term suggests that the Fermi "wound" should heal 
more quickly in the full SU(2) C calculation than in the painted case, even for large L. This 
qualitative difference in the Ansatze's structure will be manifest in a comparison of the 
models' two-body densities. 



V. MONTE CARLO COMPUTATIONS 

In this section, we shall discuss the techniques necessary to evaluate the integrals in 
Eq. ( |3.10| ) and Eq. ( |3.12| ), so that the total energy of the system and its two-body density 
may be computed. First, we shall discuss the computation of the system's energy. Let us 
consider, then, the evaluation of the kinetic energy: T = — (1/2) J2i df- Both Eq. ( |3.1(J| ) and 
Eq. ( |3.12|) have terms purely of the form 



/r2345678.../l2345678... > (5-1) 

where 

A A 2345678... = exp[-A[(xi - x 2 ) 2 + (x 3 - x 4 ) 2 + (x 5 - x & ) 2 + (x 7 - x 8 ) 2 H ]] (5.2) 

and /i2345678... nas been introduced in Eq. ( |3.5|) . For wavefunctions of this structure, the 
kinetic energy can be calculated in a straightforward way. Note that 
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(5.3a) 



whereas 

3l/l A 234... = -<9 2 /l A 234... , (5.3b) 

so that Yh=i difi2M..Pifi234... = 0- This cancellation occurs pairwise, so that 

i JV i N 

~«E ^(/l A 234.../l^4...) = E[(^/l A 234...) A^4... + /iV.^J (5.4a) 

z i=i z i=i 

= [AiV - 4A 2 ((a; 1 - x 2 ) 2 + (a* - x 4 ) 2 + ■ ■ 0)]/£ B *.„/3L (5.4b) 

A fFG 
1234.. .i 1234... 



+t F G^ /l234.../ F< " 



Note that £fg is the Fermi gas kinetic energy per quark for a finite number of quarks. 
Upon choosing anti-periodic boundary conditions, the allowed momenta in the box are odd 
integers, and 

2tt 2 N/4 

*fg P = t^E(2*-1) 2 , (5-5) 

i=l 

so that 

vrV 4 

47 = (1 - jp) . (5-6) 

where the density p = N/L has been introduced. For large N, this is the continuum SU(2) 
Fermi gas kinetic energy; the corrections to this limit are of order 1/N 2 . We hope that this 
effect will persist for other observables and for finite A. 

With the use of Eq. (|5.4b ) and Eq. (|5.6|) , the form of the derivative terms can be realized 
analytically. We are able, then, to write down the precise form of the matrix elements we 
wish to calculate. Presuming that the initial matrix element is of the form given in Eq. (ft.!2|) 
and letting A = (tffi,) 

TIU-> tne kinetic energy per particle (t)\ and the potential energy 
per particle (v)\ become 



and 



where 



(vh = - /Kr;Mr 4 , (5.8) 
V 1 N J{d Xi }A 2 ' V ; 



[(Xi - X 2 f + (x 3 - X 4 ) 2 + (X 5 - X 6 ) 2 + (X 7 - X 8 ) 2 H ]/l2345678... 

-[(Xi - X 4 ) 2 + (X 3 - X 2 f + {X 5 - Xq) 2 + (X 7 - X 8 f H ]/l4325678... 

+ ((N/2)\ -2) terms (5.9) 
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and V denotes the usual many-body potential, Eq. (|2.3|) . Note that the total energy per 
particle is (e)\ = (t)\+ (v)\. The numerator of Eq. Q5.7D no longer has a A 2 structure; 
indeed we "measure" g/A with respect to the distribution A 2 . Note that configurations 
with A = are generated with zero weight, so that the procedure we effect actually works 
rather well. 

Before we proceed, we should mention how the potential, Eq. (|2.3| ), is calculated. The 
economics literature [TJJ abounds with numerical algorithms for minimization problems - 
we need only recall the "travelling salesman" problem. For our purposes, determining the 
optimum pairing in one spatial dimension is simple. The continuous system is modelled by 
quarks on a ring; a particular quark on that ring can be paired with its neighbor to its left 
or its right, and that one choice determines the partners for all the other quarks on the ring. 
Thus, we only have to calculate the energy of the two different arrangements to determine 
the optimum pairing. Clearly, any arrangement which involves overlapping strings can not 
be an optimum pairing. 

Our results depend on the variational parameter A, and we wish to determine the optimal 
A in an efficient way. We must find X t such that 



dA 



. 



(5.10) 



A=A t 



We can do this by exploiting how the matrix elements scale in A and L. The wavefunction 
Ansatz itself is a function of the dimensionless variable XL 2 alone, but the final matrix 
elements do have a manifest L dependence. That is, 



(v)x ~ 



r 1 



L 2 



(5.11a) 
(5.11b) 



The t and v matrix elements' dependence on L differ, so that the optimal L - and, hence, 
the optimal A - can be computed directly from d(t)\/dX and d(v)\/d\. In this way, a 
single calculation for fixed XL 2 suffices to determine the optimal variational parameter and 
the observables for that A. We calculate d(t)\/dX and d(v)\/dX directly, as computing 
the numerical derivative in A of the stochastically calculated t and v matrix elements is 
unwieldly. Moreover, sampling the appropriate derivatives directly does not require much 
additional effort. Thus, we compute 



d 

dX 



1 



8A^ 
N 



8A / {dxi}Ag AX 2 J {dxi}(Ah + g 
iV J{dxi}A 2 + ~N J{dxi}A 2 

f{dxj}Ag 
lidx^A 2 



7,2 \ 



(5.12) 



and 



2 J{dxi}AVg 2 J{d Xi }Ag J{dxi}AVA 



N J{d Xl }A 2 N 



Kdx^A 2 ] 2 



(5.13) 
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where g is given by Eq. ( |5.9|) and 



Jl = [(Xt - X 2 f + (X 3 - X 4 ) 2 + (x 5 - X 6 ) 2 + (X 7 - X 8 ) 2 H ] 2 /l2345678... 

-[(Xi - X 4 ) 2 + (x 3 - X 2 ) 2 + {X 5 - Xq) 2 + (X 7 - X 8 ) 2 H ] 2 /l4325678.. 

+ ((iV/2)! - 2) terms . 



(5.14) 



As in Eq. (|5.7|), we measure g/A, as well as h/A and (g/A) 2 , with respect to the distribution 
A 2 . We can now proceed to construct the optimal variational solution for some fixed XL 2 . 
We shall label the observables computed at some p an d A by "0" , and the observables at 
the desired scaled point by "f . That is, we want p t and A t such that 



d(t) t d(v) t 







so that 



d(t) 
dX 



+ 



Po 
Pt 



1 1 



d(v)o 
dX 



. 



Thus, the optimal L - or, rather, p t - for fixed XL 2 is given by 



Pt = Po 



d{t)o 
d\ 

d(v)o 
d\ 



-1/4 



Introducing 



(5.15) 



(5.16) 



(5.17) 



we can compute the optimum A: 



Po 
Pt 



Xr 



A t = - 
V 



as well as the kinetic and potential energies at the new point: 

(*>o 



(*>* 



V 



(5.18) 



(5.19) 



(5.20) 



and 



(v)t = v(v)o ■ (5.21) 

Of course, (e) t = (t) t + (v) t - These scaling relations serve to generate the optimal variational 
solution for some fixed AL 2 . Their utility depends on the energy having a smooth, parabolic 
behavior in the region of X t , pt- Even if this constraint is satisfied, as it is here, we can not 
strictly say that we have determined the best variational solution for a fixed density. For 
a fixed density, there exist two separate parametric families - those with finite A and those 
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with A = 0. Scaling relations also exist for the A = family. The A = solution may have 
a lower energy than the finite A solution for some fixed density; the relative energies of the 
two families must be determined "empirically". 

In principle, we have two distinct strategies, Eqs. (|3.10| ) and (|3.12|) , for computing the 
observables associated with the Hamiltonian of Eqs. ( |2.2| ) and ( |2.3| ). The wavefunction in 
Eq. (|3.10|) is much simpler, but we must compute its normalization explicitly. We have 
compared the two schemes for both 8 and 12 quarks, and have found Eq. ( |3.12| ) to be the 
method of choice. The two strategies do yield identical results within error bars, as they 
should, but many more integration points are required to yield good results in the case 
where we must calculate the normalization. Consequently, the formulae and results we 
present exploit Eq. (|3.12|) . 

Thus far, we have described the procedure by which the optimal variational calculation 
of the energy is effected. However, we are interested in more than the energetics of the 
system; we would also like to study the structure of the ground state. We are particularly 
interested in the correlations induced in the ground state as a result of the internal color 
degree of freedom and global antisymmetry. Thus, we shall now describe how the two-body 
density may be calculated as well. The two-body density operator, p2(r), is defined as 

P2{r) = j'£ l 6{r-{x i -x j )). (5.22) 

This operator is manifestly symmetric, so that we may compute the resulting matrix element 
in the manner indicated by Eq. ( |3.12j ). We choose to normalize the matrix element of the 



two-body density by the square of the one-body density; that is, we shall compute the matrix 
element of p2{r)j p\. With this choice of normalization, the low-density limit of this matrix 
element is 

(p 2 (r)) Po ->o = |^isoi(r)| 2 < 5 
Po Po 

l^isoiMI 2 i s t ne square of the isolated hadron wavefunction; here 

/ i \i/a 
V2tt 

so that the low density limit of the clustering parameter A is 



* isol (r)| 2 = I^H e-^ r \ (5.24) 



A(p^0) = -^~.35. (5.25) 



The matrix element of Eq. (|5.22|) measures the correlation of any two quarks as a function of 



their separation r in the system's ground state. At low densities, the quark-quark correlation 
function is controlled by the isolated hadron wavefunction, as shown in Eq. (|5.23| ). At high 
densities, the system becomes a non-relativistic Fermi gas. In this limit, the two-body 
density is 

= ! _ '| M , (5.26) 
Po z 
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where, in SU(2) C , the Fermi momentum hp is given by kp = npo/2. There is a 50% 
probability that any two quarks in SU(2) C will have different color states; hence, the two- 
body density at r = is 1/2. 

The Metropolis Monte Carlo procedure we use allows us to calculate the two-body density 
in a straightforward way. That is, we generate the quarks' coordinates according to the 
distribution K^triaOTlTi-l 2 ; m or der to calculate matrix elements according to Eq. ( |3.12|) . 
We can compute the matrix element of Eq. ( |5 . 2 2| ) with r by simply tabulating the occurrences 
of X ^ CC j 111 the generated distribution over the interval [0, L] in bins of Ar. The resulting 
histogram is a discrete representation of the two-body density with r. As in the case with 
the energy, scaling can be used to calculate the two-body density for the optimum value 
of the variational parameter at some fixed density in an economical manner. The quantity 
(p2(^))p -»o/Po i s hself L- independent, but the scale in r over which it is plotted does depend 
on L. The step-size in r upon scaling from p to p t transforms according to 

(Ar), = (Ar) . (5.27) 

In this manner, the two-body density at the scaled values of the density and variational 
parameter can be computed from some initial po and Aq- 



VI. RESULTS 

In this section, we shall present our variational Monte Carlo results, computed using 
Eq. (|3.12|) . We have computed the energy and correlation function as a function of density 



for eight and twelve quarks. The low and high density limits of the model are known exactly, 
so that it is useful to review those results before proceeding to a discussion of the simulations. 
At low densities, the energy per quark e and the variational parameter A are given by the 
isolated hadron values. The isolated hadron wavefunction in our units is given by Eq. Q5.24 ), 
so that at low densities, 

(c)x - ^= - -35 (6.1a) 
(t)x - (v)x - - - 18 • ( 6 - lb ) 



Equation ( |6.1b| ) is a consequence of the virial theorem in the isolated hadron limit. The 



isolated hadron variational parameter is given in Eq. (|5.25 ). 



We shall first examine the energy and variational parameter as a function of density 
before turning to the correlation function results. 



A. ENERGY VS. DENSITY 

Figure shows the variational parameter A as a function of density for eight quarks. 
The indicated points are the lowest energy solutions as a function of density for the finite 
A parametric family. At low density, the scaling solution approaches the isolated hadron 
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limit value, Eq. ( |5.25|) . At high density, the energies of the optimal solutions of the finite A 
and zero A families cross; the variational parameter drops to zero for p ~ 1. The manner in 
which this level crossing occurs is shown in Fig. |3]. At low density, the energy as a function of 
the variational parameter has a single extremum. As the density increases, a double-welled 
structure develops; at sufficiently high density, the zero A solution is lower in energy The 
two extrema are well-separated in A; this generates the sudden drop in A with p seen in 
Fig. £| This behavior is a non-trivial dynamical consequence of the model. 

The meaning of "high" density is not yet clear, as the scale in Fig. |2| is in dimensionless 
units. In the units of this paper, the rms separation of the quarks in an isolated ground 
state hadron (noting Eqs. ( |5.24[ ) and ( |5.25| )) is 



(r 2 > 5 = - 7 =^.84 (6.2) 

^isol 

so that one hadron in an isolated hadron length implies a quark number density of 2.4. The 
number density of realistic nuclear matter at saturation is .17 fm -3 . The charge radius of 



the proton in the non-relativistic quark model |20| is .5 fm, so that putting one proton in 
the volume given by that radius corresponds to a number density of 1.9 fm~ 3 . Consequently, 
the density .95 in our units corresponds roughly to 4 times nuclear matter density, although 
one must be cautious in assigning dimensions in this simple model. 



FIG. [2|. The variational parameter A as a function of p for eight quarks. The points indicate 
the optimum finite A solutions as a function of density, and the associated errors bars are 
indicated. The hatched region indicates a transition region from finite to zero A, and it 
width is 2<j, where sigma is the error assigned to the transition density estimate. 
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FIG. The energy per particle e(A) as a function of A for various, fixed p. The lines are 
merely to guide the eye. The solid lines are e(A) versus A in the transition region to zero A; 
the densities for those curves from bottom to top are .90, .95, and 1.0, respectively. 

Figure |] shows the energy per particle of the eight quark system for the optimum vari- 
ational parameters of Fig. |2]. The energies of the finite A and zero A solutions are plotted 
in the same figure. As previously illustrated in Fig. |3], the energies of the finite and zero 
A solutions cross at p ~ 1. As the density increases, the potential energy goes to zero. In 
that limit, the total energy per particle is that of a free Fermi gas. This is natural - the 
interquark separation, and, hence, the string length, approaches zero as p — > oo. At low 
densities, the total energy per particle approaches the isolated hadron value, and the virial 
theorem (t = v) appears to be satisfied in that limit as well. Note that nuclear matter in this 
eight quark calculation, at least, is not bound. At moderate densities, a dynamical interplay 
exists between the potential and kinetic energies. That is, for a slightly larger average string 
length, the shape of the spatial wavefunction can "relax" and reduce its kinetic energy. 

We shall now consider the variational parameter A and the energy per particle as a 
function of density for twelve quarks, as we would like to see how sensitive our results are 
to finite size effects. Figure |5| shows the twelve quark results for the optimal variational 
parameter A as a function of p. The eight quark results from Fig. |2| are plotted as well for 
comparison. The twelve quark calculation approaches the isolated hadron limit as p — > 0, 
and it undergoes a transition to zero A in the same region of p. However, the results 
at intermediate density are different. The optimum variational parameter in the twelve 
quark case decreases more quickly as the system moves away from the isolated hadron limit, 
although it does not have the dip at p ~ .75 that exists in the eight quark case. If we now 
turn to the energy per particle as a function of density, we see the finite size effects are 
substantial there as well. These results are shown in Fig. |6|. The energy per particle for the 
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FIG. The total energy per particle e(p) as a function of p for eight quarks. The kinetic 
energy per particle is denoted by "f , whereas the potential energy per particle is denoted 
by "f " . The lines are merely to guide the eye. The points connected by lines are the energies 
associated with the finite A variational solution. The energies e, t, and v associated with 
the zero A solution are plotted as well; the symbol definitions are indicated in the figure. 
The error bars associated with these points are plotted, when they are large enough to be 
visible. The t(X = 0) points correspond to a free Fermi gas, so that there are no associated 
errors. 

twelve quark case obeys the high and low density limits of the model: the system evolves 
to the isolated hadron energy per particle and to equal t and v as the density goes to zero, 
and it approaches the Fermi gas limit as the density goes to infinity. Figure |5] also shows 
the level crossing behavior of the finite and zero A families seen in Fig. § and inferred from 
Fig. |5]. The "softening" of A seen in the twelve quark case at p ~ .3 is reflected in the 
slower increase in t in that density regime. The consequence is that the energy per particle 
skims along at the isolated hadron energy per particle for a non-trivial interval in p before 
finally increasing. This is in contrast to the eight quark case, where e increases steadily as a 
function of density. Neither the eight nor twelve quark case supports a bound state; yet, the 
change in e in going from eight to twelve quarks in the p ~ .3 regime suggests that a bound 
state can not be ruled out in this model as the number of quarks goes to infinity. The finite 
size effects are still larger in the p ~ .75 regime. There, the interplay between t and v seen 
in the eight quark case is simply absent. This is entirely consistent with the behavior of A 
in this density regime. The twelve quark A decreases smoothly to the transition region - the 
dip behavior seen in the eight quark case is absent. 
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FIG. [5|. The variational parameter A as a function of p for twelve quarks. The 
eight quark results from Fig. Q are shown as well, so that the the finite size effects 
may be seen more clearly. The twelve quark data is denoted by "x". The points 
indicate the optimum finite A solutions as a function of density, and the associated 
errors bars are indicated. The hatched region indicates a transition region from 
finite to zero A, and it width is 2a, where sigma is the error assigned to the tran- 
sition density estimate. This estimate is identical to that of the eight quark case. 



FIG. |6|. The energy per particle e(p) as a function of p for twelve quarks. The 
notation is as per Fig. || 
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FIG. |7[ The optimum variational parameter A in the painted model, 
Eq. ( |4.3|) , as a function of density for eight and twelve quarks. The 
curves are merely to guide the eye. The dashed curve connects the eight 
quark calculations, whereas the solid curve connects the twelve quark ones. 



FIG. |5| The energy per particle as a function of density in the painted model 
for eight and twelve quarks. The curves are merely to guide the eye. As in Fig. [j], 
the dashed line connects the eight quark calculations, and the solid line connects the 
twelve quark ones. The statistical errors that are not shown are negligible relative to 
the scale of the figure. The remaining notation is as in Fig. |j. 
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We can now compare the results of the full SU(2) calculation, which we have 
just presented, to the results of the approximate SU(2) C treatment - the "painted 
model" described in Sec. TV]. Figure [?| shows the optimum variational parameter in the 
painted model as a function of density for eight and twelve quarks. Here, as above, the 
variational parameter approaches the isolated hadron limit at low densities. At high 
densities, however, A flows smoothly to zero; there is no abrupt transition to zero A 
in this model. The finite size effects are small for they only become significant 

at higher densities. If finite size effects are to exist in an isolated density regime, 
then it is sensible that they should exist only at larger densities. The variational 
parameters of the models do not bear closer comparison; the Ansatze from whence 
they are derived are decidedly different. 

The energy per particle in the painted model is shown as a function of density in 
Fig.| for eight and twelve quarks. The model results approach the isolated hadron 
limit at low density and the free Fermi gas limit at high density. As above, the 
potential energy tends slowly to zero as the density increases. Here the finite size 
effects are small - even at larger densities. The total energies per particle in the 
magic and painted models are compared in Fig.|9|. The twelve quark calculations in 
the two models are fairly similar, even though the physics of the models is rather 
different. In comparing Figs.[| and |8|, we note that differences in t and v exist, but 



FIG. The energy per particle in the magic (full SU(2)) and painted (approximate 
SU(2)) models as a function of density for eight and twelve quarks. The starred data 
points joined by the solid curve are the twelve quark painted model calculations from 
Fig. p. The other solid curve connects the twelve quark magic model calculations 
from Fig. || The dashed curve joins the corresponding eight quark calculations from 
Fig. H. Note that the energy per particle in the full SU(2) calculation is for the 
optimal variational parameter solution - A is not required to be finite. 
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they tend to compensate. The similarity of the total energy per particle in the twelve 
quark cases does tend to support the use of the painted model for calculations of this 
ilk. It would be interesting to see if this similarity survives in calculations with more 
quarks. 



B. GROUND STATE STRUCTURE 

We shall now turn to a study of the system's ground state structure as a function 



of density. Figures |10| and |TT] illustrate the two-body density versus the interquark 
separation r for various one-body densities in the eight and twelve quark cases. Both 
the calculations in the full and approximate SU(2) C models are shown; we hope that 
their comparison will yield further insight into the structure of the ground state. 



FIG. |TD| . The two-body density as a function of the interquark separation r for 
a variety of densities for the eight quark case. Both the results of the full (magic) 
and approximate SU{2) (painted) calculations have been plotted. The full SU (2) 
points are joined by a solid line to guide the eye. The dashed curve is the two-body 
density appropriate to a free Fermi gas (Eq. (|5.26|) ), and the dotted curve shows that 
appropriate to the isolated hadron limit (Eq. ( |5.23| )). Note that in (a) p = .25, in (b) 
p = .50, in (c) p = .95 and A is the optimum finite solution, and in (d) p = .95 and 
A = 0. 
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FIG. [IT]. The two-body density as a function of the interquark separation r for a 
variety of densities for the twelve quark case. The notation is identical to that of 
Fig. UJ. 

We shall consider the eight quark results, Fig. |TD|, first. Figure [II] shows the two- 
body density from low to high ((a)— >(d)) density. Both the magic and painted models 
approach the isolated hadron limit at low density, as they should, although the painted 
model approaches the limit more quickly as p decreases. Note that in all cases P2{f) 
approaches one for sufficiently large r. The two-body density for the full SU(2) C 
case is shown at still lower density in Fig. |T2](a) ; pi{r) does approach the isolated 
hadron limit at low density, as claimed. The painted model's more rapid approach to 
the isolated hadron limit can be readily understood through analyzing the models' 
Ansdtze. Equations (|4.8| ) and (|4.9| ) show the structure of the two ground states as 



two quarks of the same color approach each other. The healing of the Fermi wound 
in the magic model contains an additional //-independent term; this is likely why 
the painted model approaches the isolated hadron limit more quickly. At moderate 
density, p = .5, the full and approximate SU(2) calculations are quite different, 
although there is no reason why they should be the same. At still higher densities, the 
two models are again similar, but that is because they are approaching the Fermi gas 
limit. At p = .95, A in the painted model is small, but finite, and the corresponding 
two-body density is nearly that of a free Fermi gas. The same is true of the magic 
model, although A exhibits a rapid drop to zero in this density regime. The finite 
A two-body density in (c) is only slightly different from the zero A result in (d); 
apparently the structure of the state changes little across the "transition." The 
result in (d) is essentially identical to the free Fermi gas result; this implies that 
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the finite-size effects are under some measure of control. We can explore this issue 



directly by studying the twelve quark results in Fig. [1 1 



The finite-size effects in P2( r ) under the change from eight to twelve quarks are 
significant in the magic model, although they are modest in the painted model. Note 
that the zero A two-body density in (d) shows reasonable agreement with the free 
Fermi gas result: the finite-size effects are in the finite A calculation. Perhaps the 
Ansatz of Eq. ( |3.6D is ineffective in describing the additional correlation the quarks 
feel due to confinement at moderate density. In particular, the two-body density now 
changes as A changes from its optimum finite value to zero at the transition density. 
In general, p2{r) in the magic model has additional strength at r = 0, compared to 
the eight quark results. This calculation still appears to satisfy the isolated hadron 
limit, however, as shown in Fig. [T^(b) . The ground state structure of the painted 
and magic models at the twelve quark level appears to be rather different. As in the 
eight quark case, the painted model approaches the isolated hadron limit much more 
rapidly. The differing structure of the two models' Ansatze, Eqs. (|4.8|) and Q4.9|), as 



FIG. [T^. The two-body density as a function of the interquark separation r for the 
full SU(2) C calculation - the "magic" model - in the low density limit. Case (a) is for 
eight quarks with p = .08, and case (b) is for twelve quarks with p = .1. In each case, 
the dotted line corresponds to the appropriate isolated hadron two-body density. 

two quarks of the same color approach each other, is likely responsible for this effect. 
The nature of the correlations included in the painted and magic model are differ- 
ent, and the two-body density is apparently quite sensitive to those differences. It 
would be interesting to see whether this sensitivity would be preserved in calculations 
performed with still more quarks. 

The approach of the magic model to the isolated hadron limit, as shown in Fig. |T2] , 
serves to illustrate another point concerning the structure of the state. At small r, the 
Monte Carlo calculations of the two-body density track the isolated hadron results 
fairly closely. At larger r, however, there is a significant "filling- in" of the two-body 
density relative to the isolated hadron result. This is in contrast to the painted model 
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calculations at low density, as depicted in Figs. |10|(a) and 0(a). There the two-body 
density tracks the isolated hadron result until it reaches essentially zero; it eventually 
picks up again because other hadron clusters exist at finite density. The difference 
in the behavior of the two models as r increases in the low density limit arises from 
the exchange terms included in the magic model. For any fixed, finite density, there 
exists some r for which the exchange term to recouple a quark pair into different 
hadrons is finite. Such terms do not exist in the painted model as the colors of the 
quarks are fixed. This effect explains the "wings" seens on the two-body density in 
the magic model shown in Fig. [T2|. 

We can understand the difference between the magic and painted models at low 
densities in yet another way. That is, the large dip below one in the painted model's 
P2{t) at large r can be interpreted in terms of a repulsive effective force between 
hadron clusters in this model. This is simply due to antisymmetry. A red quark in 
one cluster must be anticorrelated with the red quark it sees in an adjacent cluster. In 
the magic model, the inclusion of exchange terms allows some additional attraction 
to exist, since the colors of the quarks can change. As a red quark in one cluster 
approaches a red quark in an adjacent cluster, the first quark can flip its color and 
pair off with the red quark it approaches. In this way the effective attraction - relative 
to the painted model - is generated. 

We have now completed our survey of results in the full SU{2) C model. 



VII. SUMMARY AND OUTLOOK 

Here we have considered nuclear matter in the context of a non-relativistic con- 
stituent quark model with an internal color degree of freedom. In part, we have been 
motivated by a desire to understand traditional issues, such as nuclear matter binding 
and saturation, from a quark model viewpoint. Our primary interest, however, is in 
the properties of the system as a function of density. The model system we consider 
becomes a gas of isolated hadrons at low density, and a non-relativistic free Fermi 
gas at high density. Thus, it provides a natural laboratory in which to understand 
how the combined effects of confinement and antisymmetry determine the structure 
of the system as it evolves between the two limits. In particular, we wish to gain 
insight into a possible "deconfinement" phase transition at high density. 

The model we consider is a simple one. It is derived from the combined appli- 
cation of the adiabatic and strong coupling limits to QCD, and it consists of the 
following prescription: the quarks are always paired to the lowest energy set of flux 
tubes. Consequently, the model contains no dynamical gluon degrees of freedom; the 
dynamics are those of quark exchange. We remind the reader that the global min- 
imization prescription which defines our model also guarantees that it is free from 
residual color Van der Waals forces - an advantage over two-body potential models 



nj. The model is non-relativistic, although we are, in principle, interested in density 
regimes in which relativistic effects could be important. We feel that the model's sim- 
plicity and intrinsic appeal outweighs its disadvantages for the qualitative purposes 
of our study. Its primary advantage is that it places quark and hadron degrees of 
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freedom on an equal footing. The structure of the system at intermediate densities is 
due entirely to the model dynamics, namely of confinement and antisymmetry, rather 
than to the heuristic approximations obliged in a more complicated model. 

The quark exchange model has been used in previous studies of nuclear matter 
P~PH|. Our contribution here has been the complete incorporation of an internal 
SU(2) color degree of freedom, although we have continued to ignore the quarks' 
spin. By "complete", we mean that we have incorporated the internal degree of 
freedom without giving the quarks fixed colors. This is important: QCD requires 
that a many-quark system be an overall color singlet, but it does not determine the 
colors of the individual quarks. We have constructed an explicitly antisymmetric 
variational Ansatz with the internal degree of freedom, and then proceeded with a 
Monte Carlo calculation to determine the properties of the system. 

We have found that the SU(2) C results look rather different from earlier "spin 
fermion" studies Indeed, the simulations support an abrupt transition from 

finite to zero variational parameter, which describes the system's clustering into 
hadrons, at moderately high density. The system with zero variational parameter 
is a free Fermi gas. The change of the structure of the system across this transition is 
sensitive to finite-size effects. At the eight quark level, the two-body density changes 
little as the variational parameter changes from its optimum finite value to zero at the 
transition density. However, the same calculation performed at the twelve quark level 
shows a significant effect. A definitive conclusion as to the structure of the system 
as it goes through the transition regime must be reserved, then, for later study. This 
abrupt change in the variational pararmeter as a function of density is also absent 
in the "painted" - the quarks have fixed colors - SU(2) study we considered here 
for comparison. It would be interesting to see if this effect persists in three spatial 
dimensions and in full SU (3) color. 

The energy per particle as a function of density in this model shows rich behavior 
as a function of density. At moderate densities, the system is capable of tolerating 
longer string lengths in order to reduce its kinetic energy. We should note, in passing, 
that in a model with linear, rather than harmonic, confinement these effects would 
likely be more pronounced. The energy per particle rises more slowly with density 
than those of earlier studies without color, but nuclear matter in this model is not 
bound. The energy per particle does become considerably flatter with density under 
the change from eight to twelve quarks, however. The change is, in fact, big enough 
that the existence of a bound state can not be ruled out as the number of quarks 
increases. Yet, the lack of a bound state is hardly a surprise. The dynamics in this 
model are simple and do not suffice to reproduce the known features - such as the 
spin dependence - of the N — N force. 

Unlike the full SU{2) model, the painted SU(2) calculation of the energy per par- 
ticle shows small finite size effects under the change from eight to twelve quarks. The 
twelve quark energy per particle results in the full and painted SU(2) calculations are 
surprisingly similar. This is in contrast to the evolution of the variational parameter 
and two-body density of the two calculations with density, which are different. It 
would be interesting to see whether this similarity persists in a calculation with more 
quarks. The similarity which does exist, however, lends support to the further use 
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of the painted model in exploring the energetics of bulk quark matter, for example, 



when flavor degrees of freedom are included pl |. However, the distinctly different 



two-body densities in the two models indicate that a complete treatment of antisym- 
metry is necessary to understanding the system's ground state structure. Indeed, this 
may also indicate that a complete incorporation of spin and flavor degrees of freedom 
as well are essential to an understanding of the ground state of hadronic matter with 
density. 

The computation of the two-body density in this model has given us insight into 
the structure of the ground state as a function of density. In principle, we would also 
like to study the excitations of the system as well. We can explore the "phonon-like" 
behavior of the system through the calculation of the Fourier transform of the two- 
body density P2] . It is possible, however, to gain access to the full response of the 
system. That is, we can calculate the quark response function in imaginary time, and 
then reconstruct the longitudinal response function, familiar from inclusive electron 
scattering, by finding the inverse Laplace transform of the Euclidean response. The 



inversion procedure is non-unique, but Carlson and Schiavilla |23[ have shown it 
practicable for light nuclei. We imagine that the quark response function would 
evince rich behavior SIS db function of density. 
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FIGURES 



FIG. 1. An illustration of quark exchange dynamics for 4 quarks. Two incoming hadrons with 
quark content (12)o and (34)o may exchange quarks - the strings connecting 1 to 2 and 3 to 4 flip 
- to yield a (13)o(24)o outgoing state. 

FIG. 2. The variational parameter A as a function of p for eight quarks. The points indicate 
the optimum finite A solutions as a function of density, and the associated errors bars are indicated. 
The hatched region indicates a transition region from finite to zero A, and it width is 2a, where 
sigma is the error assigned to the transition density estimate. 

FIG. 3. The energy per particle e(A) as a function of A for various, fixed p. The lines are 
merely to guide the eye. The solid lines are e(A) versus A in the transition region to zero A; the 
densities for those curves from bottom to top are .90, .95, and 1.0, respectively. 

FIG. 4. The total energy per particle e(p) as a function of p for eight quarks. The kinetic 
energy per particle is denoted by "t", whereas the potential energy per particle is denoted by 'V. 
The lines are merely to guide the eye. The points connected by lines are the energies associated 
with the finite A variational solution. The energies e, t, and v associated with the zero A solution 
are plotted as well; the symbol definitions are indicated in the figure. The error bars associated 
with these points are plotted, when they are large enough to be visible. The t(X = 0) points 
correspond to a free Fermi gas, so that there are no associated errors. 

FIG. 5. The variational parameter A as a function of p for twelve quarks. The eight quark 
results from Fig. || are shown as well, so that the the finite size effects may be seen more clearly. 
The twelve quark data is denoted by "x". The points indicate the optimum finite A solutions as 
a function of density, and the associated errors bars are indicated. The hatched region indicates a 
transition region from finite to zero A, and it width is 2a, where sigma is the error assigned to the 
transition density estimate. This estimate is identical to that of the eight quark case. 

FIG. 6. The energy per particle e(p) as a function of p for twelve quarks. The notation is as 
per Fig. |||. 

FIG. 7. The optimum variational parameter A in the painted model, Eq. ( |4.3| ), as a function 
of density for eight and twelve quarks. The curves are merely to guide the eye. The dashed curve 
connects the eight quark calculations, whereas the solid curve connects the twelve quark ones. 

FIG. 8. The energy per particle as a function of density in the painted model for eight and 
twelve quarks. The curves are merely to guide the eye. As in Fig. [/], the dashed line connects the 
eight quark calculations, and the solid line connects the twelve quark ones. The statistical errors 
that are not shown are negligible relative to the scale of the figure. The remaining notation is as 
in Fig. g. 
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FIG. 9. The energy per particle in the magic (full SU(2)) and painted (approximate SU(2)) 
models as a function of density for eight and twelve quarks. The starred data points joined by 
the solid curve are the twelve quark painted model calculations from Fig. |8|. The other solid 
curve connects the twelve quark magic model calculations from Fig. || The dashed curve joins the 
corresponding eight quark calculations from Fig. |4| Note that the energy per particle in the full 
SU(2) calculation is for the optimal variational parameter solution - A is not required to be finite. 

FIG. 10. The two-body density as a function of the interquark separation r for a variety of 
densities for the eight quark case. Both the results of the full (magic) and approximate SU(2) 
(painted) calculations have been plotted. The full SU(2) points are joined by a solid line to guide 
the eye. The dashed curve is the two-body density appropriate to a free Fermi gas (Eq. ( |5.26D ), 
and the dotted curve shows that appropriate to the isolated hadron limit (Eq. (|5.23[) ). Note that 
in (a) p = .25, in (b) p = .50, in (c) p = .95 and A is the optimum finite solution, and in (d) p = .95 
and A = 0. 

FIG. 11. The two-body density as a function of the interquark separation r for a variety of 
densities for the twelve quark case. The notation is identical to that of Fig. [k| 

FIG. 12. The two-body density as a function of the interquark separation r for the full SU(2) C 
calculation - the "magic" model - in the low density limit. Case (a) is for eight quarks with 
p = .08, and case (b) is for twelve quarks with p = .1. In each case, the dotted line corresponds to 
the appropriate isolated hadron two-body density. 
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